Band gap control via tuning of inversion degree in CdIn2S4 spinel 
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Based on theoretical arguments we propose a possible route for controlling the band-gap in the 
promising photovoltaic material CdIn2S4. Our ab mitio calculations show that the experimental 
degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approxi- 
mately to the equilibrium value given by the minimum of the theoretical inversion free energy at a 
typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthe- 
sis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic 
band-gap of the solid, as shown here by accurate screened hybrid functional calculations. 



The cadmium/indium thiospinel CdIn2S4 is a photo- 
sensitive semiconductor with excellent light absorption 
properties in the visible range of the spectrum, and has 
attracted considerable attention in recent years due to its 
potential applications in photocatalysis, high-efficiency 
solar cells, light-emitting diodes, and optoelectronic de- 
vices IJ-Q- A detailed understanding of the factors con- 
trolling its electronic band structure is required in order 
to optimize these applications. The band gap of pure 
CdIn2S4 is indirect and values between 2.1 eV and 2.4 
eV (between 2.5 and 2.7 eV for the direct gap) have been 
reported by different authors Q. The electronic and op- 
tical properties of the system at room temperature seem 
to depend on whether the crystal is annealed or quenched 
from the synthesis temperature, which is possibly related 
to changes in the cation distribution Q. 

In a 'normal' spinel, the 2+ cations are located in the 
tetrahedral sites and the 3-f cations occupy octahedral 
sites, while deviation from this distribution is called 'in- 
version'. For CdIn2S4 it is difficult to refine the Cd/In 
occupancies of the tetrahedral and octahedral sites from 
standard diffraction measurements because Cd^"*" and 
In'^+ are isoelectronic, but Raman experiments suggest 
that some level of inversion is present 0, [1| . We can write 
the formula unit as (Cdi_a;In2;)[In2_2:Cda;]S4, where ( ) 
represent the tetrahedral sites, [ ] the octahedral sites, 
and X is the degree of inversion (x=0 is the normal spinel 
and x=l is the fully inverse spinel). The experimental 
value of a; is 0.20 0. 

In this Letter we present theoretical results show- 
ing that (j) the experimental degree of inversion degree 
roughly corresponds to thermodynamic equilibrium at 
the formation temperatures of the thiospinel, thus sug- 
gesting that the distribution of cations could be modified 
by changing the synthesis temperature, and (ii) that even 
a moderate change in the degree of inversion leads to sig- 
nificant change in the band gap of the material, which 
indicates a promising route for tuning the light absorp- 
tion properties of the semiconductor. 
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We have performed density functional theory (DFT) 
calculations of CdIn2S4 with different inversion degrees, 
using the Vienna Ab Initio Simulation Program (VASP) 
9]. We employed both the Perdew-Burke-Ernzerhof 
(PBE) exchange-correlation functional jl^ and its mod- 
ified version for sohds (PBEsol) which improves the 
description of cell geometries and of phonon frequen- 
cies in solids, compared to the standard PBE functional 
[l^ . For the direct CdIn2S4, PBEsol gives an equilib- 
rium cell parameter (a=10.860 A) that is closer than 
the PBE result (a=11.020 A) to the experimental value 
(ae2;p=10.831 A, which can be obtained by extrapolation 
to absolute zero of the temperature dependence measured 
by Kistaiah et al. [l^). The convergence of energies with 
respect to cutoff energies, reciprocal space sampling den- 
sity and other precision parameters was checked carefully. 
The CdIn2S4 primitive cell contains two formula units, 
which allows us to consider inversion degrees x — 0, 0.5 
and 1 . For each value of x there is only one symmetrically 
different configuration of cations, so we calculate the in- 
version energy as AEconf{x) = E{x) — -E(O), where the 
subscript conf indicates that this is the configurational 
contribution only (vibrational contributions are discussed 
below). 

A quadratic dependence of the inversion energy with x 
has been empirically recognised by Kriessman and Harri- 
son [lj| and theoretically justified by O'Neill and Navrot- 
sky |15| based on the linear dependences with x of both 
the cell parameter a and the anion position parameter u. 
Therefore, our three calculated inversion energies allow 
us to interpolate for any value of x using a quadratic func- 
tion, as shown in Fig [1] a) . The PBE and PBEsol results 
are very similar in terms of energy, therefore the discus- 
sion below refers to the PBEsol results, unless otherwise 
stated. The inversion energies are positive for the whole 
range of x, and the curvature is slightly negative. In or- 
der to test the approximation of using only one configu- 
ration for the calculation of the energy for each a; > 0, we 
have calculated the inversion energy for x=0.25 by taking 
the symmetry-adapted ensemble average of all different 
cation configurations [ll] in a cell doubled along one axis. 
The result, represented as an empty circle in Fig. [T]a), 
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FIG. 1. a) Inversion energies as obtained from PBE and 
PBEsol calculations; the open circle corresponds to the en- 
semble average in a double cell with a;=0.25. b) Configura- 
tional and c) vibrational contributions to the inversion free 
energy, d) Total inversion free energy; the vertical line marks 
the experimental value of the inversion degree. 



is in good agreement with our quadratic interpolation 
based on the primitive cell results. 

We can now estimate the configurational free energy 
of inversion AFconf = ^E^onf ~ TASconf, where: 



Sconf = —R[xlnx + (1 — x)ln{l — x) + 

xln^ + {2 - x)ln{l ^ ^)] (1) 

is the ideal configurational entropy of inversion [l7l . [l8l |. 
In practice, some excess (non-ideal) contributions to the 
configurational entropy can be expected, but our test cal- 
culations for the ensemble of configurations with x—0.25 
in the double cell show that these excess contributions are 
indeed small: the difference between the temperature- 
dependent entropy (as calculated using Boltzmann statis- 
tics, e.g. [l^) and the maximum entropy for the given 
cell and composition, is only 2%. In the absence of in- 
version energies, the minimum of the inversion free en- 
ergy corresponds to the maximum of the configurational 
entropy, which occurs at x = 2/3 (full disorder of the 
cations among all sites). However, the positive and rel- 
atively high inversion energies imply that the minima of 
the inversion free energy actually occur at much lower 
values of x (e.g. x=0.04 at T=1000 K), as shown in Fig 
[T]b). Therefore, according to the analysis so far, which 
ignores vibrational effects, the CdIn2S4 spinel would be 
expected to be almost fully direct, while experimental 
measurements suggest inversion degrees of ~20% |8j]. 

In order to discuss the effect of vibrations in the ther- 
modynamics of inversion, we now examine the phonon 
modes in CdIn2S4 with different degrees of inversion. For 



these calculations, we employ density functional pertur- 
bation theory, as implemented in VASP. We first com- 
pare the calculated zone-centre frequencies for the di- 
rect spinel with the experimental values measured us- 
ing infrared and Raman spectroscopy 0, H^. TABLE 
H] shows that the vibrational modes are well described 
by our PBEsol calculations, with average discrepancy of 
only 2.5%. 

We then calculated the phonon frequencies in a 2 x 2 x 2 
supercell, which is equivalent to a F-centered 2x2x2 
sampling of the reciprocal space, thus allowing (linear) 
dispersion of the modes. From the resulting frequencies 
we obtained the vibrational contributions to the inversion 
free energy in the harmonic approximation: 
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for configurations with inversion degrees x—0, 0.5 and 
1. Fig. [T] c) shows that the vibrational contribution to 
the free energy exhibits an almost linear variation with 
x. The negative slope means that vibrational effects 
will shift the equilibrium inversion degree towards values 
higher than those expected based on configurational con- 
tributions only. The total (configurational -I- vibrational) 
free energy is then plotted in Fig. [TJl) as a function of x 
for temperatures between 900 and 1200 K. At a tempera- 
ture T= 1200 K, which is typical for the synthesis of this 
type of material [H, [2l| , the equilibrium degree of inver- 
sion is very close to the experimentally observed value 
x=0.20. This result suggests that the degree of inversion 
in CdIn2S4 is thermodynamically controlled during the 
solid formation. 

The sensitivity of the equilibrium inversion to the tem- 
perature indicates that it is possible to tune the cation 
distribution via temperature control during the sample 
preparation. For example, if the equilibration tempera- 
ture is reduced from 1200 K to 1100 K, the degree of 
inversion decreases to x=0.15. From these results, it 
can also be expected that the degree of inversion will 



TABLE I. Calculated zone-centre phonon frequencies of nor- 
mal CdIn2S4 spinel in comparison with experimental values 



Mode 






Tin (IR) 


68 


67 


T2g (Raman) 


93 


92 


Ti„ (IR) 


171 


171 


Eg (Raman) 


185 


189 


ri„ (IR) 


215 


221 


T2g (Raman) 


247 


239 


ri„ (IR) 


307 


296 


T2g (Raman) 


312 


297 


A\g (Raman) 


366 


355 
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FIG. 2. Direct and indirect band gaps as a function of inver- 
sion degree, from screened hybrid functional calculations. In 
the inset, the band gaps as a function of the configurational 
equilibration temperature. 

be higher when the sample is quenched (cooled rapidly) 
instead of annealed (cooled slowly) after its formation. 
In the latter case, the slow cooling will allow the equili- 
bration of the cation distribution at lower temperatures. 
Experimental studies have indeed shown a variation in 
the electronic and optical properties of CdIn2S4 with the 
cooling rate after synthesis 

We therefore consider the effect that a change of inver- 
sion degree has in the electronic properties of the mate- 
rial. The electronic structure was calculated here using 
the Heyd-Scuseria-Ernzerhof (HSE06) screened hybrid 
functional, which includes 25% of Hartree-Fock exchange 



[23 |, yielding better band gap predictions (although at 
a higher computational cost) than standard DFT func- 
tionals [23]. In fact, recent work has shown that HSE06 
gives band structures for Mg/In and Cd/In thiospinels 
that are very similar to those obtained from more com- 
putationally demanding many-body techniques The 
calculated values for the direct and indirect band gap are 
close to, although somewhat lower than, the experimen- 
tal values measured in samples with degree of inversion 
a;=0.20 (2.1-2.4 for indirect and 2.5-2.7 for direct band 
gap 0). The effect of inversion on the band gap is found 
to be quite drastic, with a decrease of almost 1 eV in 
both the direct and indirect band gaps when x changes 
from to 0.5 (Fig. EJ. 

This result means that any small change in the inver- 
sion degree, achieved via temperature control during the 
sample preparation, will be reflected significantly in the 
electronic and optical properties of the material. The 
connection between the band gap and the equilibration 
temperature for the cation distribution is illustrated in 
the inset of Fig. [2] Lower temperatures of formation, or 
slower cooling rates after sample preparation, should lead 
to wider band gaps. Therefore our simulations results in- 
dicate a possible route to tune the electronic properties 
of this interestic photovoltaic material. 
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